function find_error(fname1,fname2,nx,ny)
  d1 = load(fname1);
  d2 = load(fname2);
  a1 = reshape(d1(:,3), ny, nx);
  a2 = reshape(d2(:,3), ny, nx);
  fprintf('%e\n', sum(sum(abs(a1-a2).^2)));
end